High-Mobility Few-Layer Graphene Field Effect Transistors Fabricated on Epitaxial 

Ferroelectric Gate Oxides 
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The carrier mobility /i of few-layer graphene (FLG) field-effect transistors increases ten-fold when 
the Si02 substrate is replaced by single-crystal epitaxial Pb(Zro.2Tio.s)03 (PZT). In the electron- 
only regime of the FLG, fi reaches 7xl0 4 cm 2 /Vs at 300K for n = 2.4 x 10 12 /cm 2 , 70% of the 
intrinsic limit set by longitudinal acoustic (LA) phonons; it increases to 1.4x 10 5 cm 2 /Vs at low 
temperature. The temperature-dependent resistivity p(T) reveals a clear signature of LA phonon 
scattering, yielding a deformation potential D = 7.8±0.5 eV. 

PACS numbers: 73.50.-h, 72.10.-d, 77.84.Dy 



Recent calculations show that the intrinsic mobility 
of graphene, set by longitudinal acoustic (LA) phonon 
scattering, can reach ^10 5 cm 2 /Vs at room tempera- 
ture However, extrinsic scattering sources, many 
of which arise from the surface morphology, chemistry, 
structural and electronic properties of the widely used 
SiC>2 substrate, limit the mobility to the current range 
of 2 x 10 3 - 2 x 10 4 cm 2 /Vs [H, % 0, 0, 0, 0, 0, 0, 0, El 
. Increasing the mobility beyond these extrinsic limits 
is one of the central challenges of the graphene commu- 
nity. Recently, two groups have reported a significant 
improvement in the mobility of suspended graphene af- 
ter current-heating annealing [ill . u3^ . A more device- 
friendly solution involves placing graphene on a different 
substrate. Studies on several alternatives have been re- 
ported recently although these substrates result in mo- 
bilities comparable to that on SiC>2 fl3j |. 

In this letter we report significant carrier mobility 
improvement in few-layer graphene (FLG) field effect 
transistors (FETs) fabricated with single-crystal epitax- 
ial Pb(Zr .2Tio.8)0 3 (PZT) films as the gate oxide. At 
300 K, PZT-gated FLG exhibits a mobility fi ~ 7 x 10 4 
cm 2 /Vs at a density of n = 2.4 x 10 12 /cm 2 , reaching 
70% of the intrinsic limit set by LA phonons. We ob- 
serve a clear signature of LA phonon scattering in the 
temperature dependence of resistivity p(T). The PZT- 
gated FLG shows a residual resistivity p at low temper- 
ature approximately an order of magnitude lower than 
that of Si02-gated single and few-layer graphene. This 
low po corresponds to a mobility of 1.4 x 10 5 cm 2 /Vs 
and a long carrier mean free path of 2 pm at n = 
2.4 x 10 12 /cm 2 . Our results open up a promising route 
into realizin g gr aphene's full scientific and technological 
potential pi Tld fla] . 

For our FETs we work with 400 nm PZT films that 
are epitaxially grown on Nb-doped single-crystal SrTi03 
(STO) substrates via radio- frequency magnetron sputter- 
ing [161 ] . X-ray diffraction and atomic force microscopy 
(AFM) measurements show that these films have high 
crystalline and surface quality (Fig. Ufa)) fl7j |. 



FLG flakes are mechanically exfoliated on PZT from 
Kish graphite (Toshiba Ceramics Co.) followed by optical 
identification and AFM characterization (Fig.QJa)) [l7j ]. 
Multiple electrodes (30 nm Au/3 nm Cr) in the Hall-bar 
configuration are patterned onto regularly shaped flakes 
using e-beam lithography and metal evaporation. The 
Nb-doped STO substrate serves as the backgate to which 
a bias voltage V g is applied to tune the carrier density 
of the FLG (Figs. QJb) and (c)). Results reported here 
are collected from 3 FETs fabricated on the same PZT 
substrate and one FET on a Si02 substrate. 
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FIG. 1: (a) AFM contact mode image of a 2.4 nm FLG flake 
(center) on a 400 nm PZT film. Inset: optical image of the 
whole flake with the area in (a) circled. The PZT surface 
shows smooth terraces separated by a-axis lines, with a root- 
mean-square (RMS) surface roughness of 3-4 A over a 1 /im 2 
area. FLG has a roughness of 2-3 A. (b) Device schematics, 
(c) Hall bar configuration of a FLG-FET with current (Ii, I2) 
and voltage electrodes for resistance (Vi and V2) and Hall (Vi 
and Vh) measurements. We determine the thickness of this 
FLG to be (2.4±0.3) nm based on its optical transparency. 

Resistivity and Hall measurements were performed in a 
4 He cryostat with a base temperature of 1.4 K, equipped 
with a superconducting magnet. Standard low frequency 
(47Hz) lock-in techniques are used with excitation cur- 
rents ranging from 50 to 200 nA. In Fig. [2] we show the 
sheet resistivity p of a FLG device (~2.4 nm or 7 layers, 
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Fig. [He)) as a function of V g at temperatures 4 K < T < 
300 K. p(V g ) displays a broad maximum at the charge 
neutrality point. Curves below 300 K are shifted to align 
the p(V g ) maximum at V g = V 0. FLG of this thick- 
ness behaves as a two-dimensional (2D) semimetal, where 
the low-energy bands for electrons and holes are parabolic 
and overlap slightly [19( (inset of Fig.[2ja)). The carrier 
density in the FLG is controlled by V g through n e -rih — 
aV g , where a is the charge injection rate of the backgate. 
In the band-overlap regime (regime I in Fig. [2la) inset), 
both electrons and holes contribute to conduction: 
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- = e(n e p e + n h ph) 
P 



(1) 



At sufficiently large |V S |, the system becomes a pure 2D 
electron (regime II in Fig. [D(a) inset) or hole (not shown) 
gas . There, the resistivity and the Hall coefficient Rh 
are given by: 



1 1 

- — en e ^p e ,h', Rh = 

P en e ,h 



; n e h 



aV a 



(2) 



We measure Rh in the hole-only regime of two devices 
and determine a = 1.35 x 10 12 cm _2 /V. Using a parallel- 
plate capacitor model, we extract a dielectric constant 
K ~ 100 for our PZT films. This value is confirmed by in- 
dependent low-frequency capacitance measurements (see 
Ref. [I3| for details). Compared with k = 3.9 in Si02, 
PZT gate oxides are extremely efficient in injecting carri- 
ers into graphene, as well as screening charged impurities 
in its vicinity. 
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FIG. 2: (a) p(V g ) at selected temperatures taken on the device 
shown in Fig.fjjc). Inset: schematics of the band structure of 
FLG of this thickness, (b) p(V g ) at 4 K. The kink at Vj = 1.1 
V (dash-dotted line) marks the boundary between regimes I 
and II. (c) p(V g ) at 10 K (open symbols) with a fitting curve 
(solid line) combining Eqs. [TJ and [3] with f3 — 0.9 and r = 
0.6. The dashed line is calculated from Eq. [T] assuming a 
density-independent mobility fi e — fih = 1 x 10 5 cm 2 /Vs. 



It is clear from Eqs. [T] and [2] that the slope of p(V g ) 
changes at a threshold Vj ', where the system transitions 
from a two-carrier to a single-carrier regime. Such a 



slope change has been reported in Bismuth nanowires, 
which is another two-carrier system [20(. Indeed, a kink 
at Vj = 1.1 V is clearly visible in p(V g ) at low tem- 
perature (Fig. [2£b)), where n e — 1.5 x 10 12 /cm 2 and rih 
= 0. Modeling the FLG in regime I with one electron 
and one hole band and using the effective mass values 
determined in Ref. [lj| for this thickness (m* = 0.06 
mo and m* h = 0.10 mo), we estimate the electron and 
hole densities at the charge neutrality point to be n° 
= »J ~ 9 x 10 11 /cm 2 . This corresponds to an overlap 
between the electron and hole bands of ~30 meV (see 
Ref. 17| for more discussions). These estimates are in 



good agreement with results obtained using methods de- 
scribed previously fl9| and band structure calculations 
of FLG of this thickness 21]. These studies also suggest 



that FLG in this thickness range may have more than 
one hole band 

0, l2ll ]. We emphasize that the central 
results of the present study are taken in the electron-only 
regime described by Eq. [21 and do not rely on the accu- 
rate knowledge of the band structure in the two-carrier 
or hole-only regimes. 

In single and few-layer graphene fabricated on Si02 
substrates, the mobility has been found to be roughly 
density- independent and explained in terms of long-range 
charged impurity scattering @, i, 0, i, H, 0, p(V g ) 
calculated using this assumption and also using Eq. [1] 
is plotted in Fig. 0(c) (dashed curve). Clearly the 
n-independent assumption does not describe our data 
(open symbols) in the band-overlap regime (I). Instead, 
a density-dependent mobility p e h ~ . produces an 
excellent fit to the data within the entire regime. The 
power-law functional form is motivated by measurements 
in regime II, shown later. The solid line in Fig. [2](c) shows 
such fitting with mobilities determined by: 



Pe(n e ) 



P^h(n h ) 



(3) 



where we require p e and p,h to have a power-law depen- 
dence on n e and nh respectively with the same exponent 
P but scales by a factor r. We obtain f3 = 0.9 from the 
fit in Fig. dfc). The constant c is determined by match- 
ing a measured data point /x e = 1.0 x 10 5 cm 2 /Vs at the 
electron density of n = 1.75x 10 12 /cm 2 in regime II. The 
approximate symmetric V^-dependence p(V g ) displayed 
for both carriers in regime I, together with /3 ~ 1, im- 



plies that r = 



m* e /m* h 



0.6. In addition to 



p{V g ), the above fitting parameters also produce excel- 



lent agreement with Rh data in regime I (see Ref. [17 1 
for details). The origin of the mobility asymmetry is un- 
clear to us at the moment although we note that such 
phenomenon also occurs in some Si02-gated graphene 
devices [ll|, HI}- The n-dependence of p,, on the other 
hand, distinguishes our device from Si02-gated graphene 
and FLG BE 0, EjJ • This n-dependence may originate 
from the overlapped bands, or suggest different scatter- 
ing mechanisms. Results from the electron-only regime 
shown later seem to support the second explanation. 
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FIG. 3: (Color online) p(T) at electron densities of (from top 
to bottom) n = 1.89, 2.02, 2.16, 2.30 and 2.43x 10 12 /cm 2 . 
The solid lines are fittings to Eq. FJ for T > 100 K, with 
the corrections due to a non-degenerate Fermi gas included. 
Inset: Low-T residual mobility fio(n) in a double-log plot. 
Open squares are data taken in regime II. The dashed line 
plots the fitting (Eq. [3] electrons) obtained in regime I. 
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FIG. 4: (Color online) Comparison of p(T) in various 
graphitic materials. Solid squares: PZT-gated FLG shown 
in Fig. El at n = 2.4x 10 12 /cm 2 . Open triangles: a Si0 2 - 
gated FLG of the same thickness and density [T^. Open 
circles: single-layer graphene on Si02 reported in Ref. [10]. 
Crosses: mobility of bulk graphite from Ref. [2^]. Solid line: 
LA phonon-limited mobility calculated from Eq. [4] 



The above analysis provides an approximate scenario 
of transport in the band-overlap regime of the FLG. Be- 
low we present and analyze the central results of our 
work, derived from data taken in regime (II) of the FLG 
{Vj > 1.1 V), where the FLG behaves as a single- 
carrier, one band, two-dimensional electron gas described 
by Eq. [21 Figure [3] plots p(T) extracted from data 
shown in Fig. Ufa) at five electron densities ranging from 
1.9xl0 12 /cm 2 (at V g = 1.4 V) to 2.4xl0 12 /cm 2 (at V g 
— 1.8 V), well into regime II. At a fixed n, p{T) fol- 
lows a linear T-dependence between 100 K and 300 K 
and quickly saturates to a non-zero residual value po(n) 
at lower T. This linear T-dependence, its tempera- 
ture range, and the magnitude of the resistivity change 
strongly point to scattering between electrons and LA 
phonons in graphene. This phonon mode, nearly identi- 
cal in graphene and gra phi te, has been calculated [l[ and 
experimentally studied |10j ] recently in graphene on Si02- 
However, the combination of a large po and the onset of 
another scattering mechanism at 150 K in Si02-gated 
graphene makes it difficult to extract the LA phonon 
contribution unambiguously in those systems [lfj. In 
our devices, a small po enables us to clearly observe 
the predicted linear T-dependence at T > Tbc where 
Tbg — 2nk ^ ph fa 80K is the Bloch-Gruneisen tempera- 
ture at n = 2xl0 12 /cm 2 , using a sound velocity v p h — 
2.1xl0 6 cm/s for LA phonons in graphene and /c_f = y/wn 
for the Fermi wave vector of the 2D electron gas. At 
T > Tbg , the contribution to the resistivity from LA 
phonons is given by: 



p p h(T,n) = — !(-) = - 



1 (m* e ) 2 D 2 k B T 



n 4h 3 e 2 pr, 



(4) 



u ph 



where we have modified the derivation in Ref. [lj to ac- 
count for massive electrons in FLG. D is the unscreened 
acoustic deformation potential [13] and p rn = 6.5xl0~ 7 
kg/m 2 is the areal mass density of graphene. The cor- 
rection due to a non-degenerate Fermi gas is less than a 
few percent in our density and temperature range and is 
neglected in Eq. [U 

Solid lines in Fig. [3] show fittings at different densities 
for T > 100 K, where the slopes range from 83 to 87 
mtt/K and lead to D — 7.8±0.5 eV in graphene. This 
result falls within the range of reported values in the 
literature of 1 - 30 eV [13, 0, 0, 0, [13] and agrees 
very well with tight-binding calculations producing D ~ 
37 , where 7^3 eV is the nearest-neighbor hopping 
matrix [251 ] . We do not observe evidence of super-linear 
T-dependences reported in graphene on Si02 @, [l(| that 
are attributed to remote substrate [3, [13] or inter-ripple 
flexural phonons [7] . We speculate that a higher stiffness 
and a larger average carrier-substrate separation in FLG 
may be responsible for suppressing scatterings from these 
two types of phonons. 

The small residual resistivity po in PZT-gated FLG 
leads to mobility po in excess of 1 x 10 5 cm 2 /Vs at low 
T. Since both FLG and single layer graphene are sub- 
ject to similar scattering mechanisms, a comparison be- 
tween p of PZT-gated FLG, Si0 2 -gated FLG and Si0 2 - 
gated graphene highlights the important role played by 
the substrate. Such comparison is shown in Fig.[4l where 
we compare p{T) obtained from two 2.4 nm thick FLG 
(one on PZT, one on Si02[I3])) graphene on SiC"2 from 
Ref. [13], bulk graphite from Ref. [3] and the intrinsic 
LA phonon-limited mobility calculated from Eq. [51 us- 
ing D — 8 eV. At a density of n = 2.4x 10 12 /cm 2 , the 



4 



PZT-gated device shows /i ~ 7 x 10 4 cm 2 /Vs at room 
temperature, 70% of the intrinsic phonon mobility of 
~ 1 X 10 5 cm 2 /Vs. At low T, fi increases to 1.4 x 10 5 
cm 2 /Vs, corresponding to a long mean free path of 2 /j,m. 
A second device (~5nm thick, not shown) on the same 
PZT substrate exhibits mobilities of 7.5 x 10 4 cm 2 /Vs at 
room temperature and 1.5 x 10 5 cm 2 /Vs at low temper- 
ature. These values represent an approximately ten-fold 
increase over those of our Si02-gated FLG as well as 
single and few-layer graphene reported in the literature, 
where /i ranges 2 x 10 3 — 2 x 10 4 cm 2 / Vs with weak or no 
temperature dependence @, [7 , 13, HI ■ This remarkable 
mobility improvement clearly demonstrates the advan- 
tage of the PZT substrate over SiC>2 towards fabricating 
graphene-based high-quality 2D systems. 

The low-temperature residual mobility no{n) in PZT- 
gated FLG exhibits a density dependence best described 
by no(n) ~ n 13 for 1.9x 10 12 /cm 2 < n< 2.4 xl0 12 /cm 2 . 
In the inset of Fig. [3l we show no(n) data in this range to- 
gether with the fitting obtained in regime I: /J>o(n) ~ n 9 . 
This n-dependence of /z is in contrast to the Si02-gated 
graphene, where the scattering due to Coulomb impuri- 
ties leads to a very weak n-dependence in a comparable 
density range, suggesting that different scattering mech- 
anisms are at work @, [J i, |, 1,1,01,0. 

It has been shown in suspended graphene that a sig- 
nificant improvement in fi is only achieved after current 
annealing, which highlights the important role played by 
interfacial adsorbates [11[, among otherproposed scat- 
tering mechanisms at low temperature 

0, ill Hill. 

Our PZT substrates possess a large spontaneous polariza- 
tion P pointing into the surface [LJj . The absence of free 
carriers in ungated FLG devices indicates that this polar- 
ization is almost completely screened by a high-density 
layer of surface adsorbates prior to exfoliation. Screening 
adsorbates may come from free ions, atoms and molecules 
in the ambient and OH~ and H + produced by the disso- 
ciation of H2O on PZT surface 3l|, 32, 33|. Despite their 
high density, our data suggest that the scattering from 
interfacial adsorbates is much weaker than in Si02-gated 
devices. We attribute this remarkable phenomenon to 
the strong screening of PZT and speculate that some de- 
gree of ordering in the adsorbate layer may also play a 
role in reducing the scattering. 

In conclusion, we have demonstrated a significant per- 
formance improvement in few-layer graphene FETs by 
using the crystalline ferroelectric gate oxide PZT. This 
approach has led us to the observation of the highest re- 
ported mobility to date in unsuspended single and few- 
layer graphene devices. This result opens up a new route 
for realizing high-speed electronic devices and exploring 
novel 2D physics in graphene. 
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